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ABSTRACT 

We examine the brightness of the Cosmic Radio Background (CRB) by comparing 
the contribution from individual source counts to absolute measurements. We use a 
compilation of radio counts to estimate the contribution of detected sources to the CRB 
in several different frequency bands. Using a Monte Carlo Markov Chain technique, we 
estimate the brightness values and uncertainties, paying attention to various sources of 
systematic error. At = 150 MHz, 325 MHz, 408 MHz, 610 MHz, 1.4 GHz, 4.8 GHz, and 
8.4 GHz our calculated contributions to the background sky temperature are 18, 2.8, 
1.6, 0.71, 0.11, 0.0032, 0.0059 K, respectively. We then compare our results to absolute 
measurements from the ARCADE 2 experiment. If the ARCADE 2 measurements are 
correct and come from sources, then there must be an additional population of radio 
galaxies, fainter than where current data are probing. More specifically, the Euclidean- 
normalized counts at 1.4 GHz have to have an additional bump below about 10 fiJy. 
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1 INTRODUCTION 

Investigating what sources make up the diffuse extragalac- 
tic background over a wide range of wavelengths can help 
us to understand the different physical mechanisms which 
govern t he generation and transport of energy over cosmi c 
time (e.g iLongair fc Sunvaev|[l969l : iRessell fc Turneijll990D . 
Much effort has gone into resolving the sources which com- 
prise the backgr ound at 7-ray, X-ray, optical, and infra-red 
wave l engths (e.g.lMadau fc Pozzettil.l2000l;lHauser &: Dwekl . 



l200ll : lBrandt fc Hasingeil . l2005l : lLagache et al.1 . 120051 ). Uaw- 
ever, the radio part of the spectrum has received far less 
attention. While there have been many radio surveys and 
compilations of source counts done over the years, there 
have been only a few attempts at using t hese to obtain 
estimates of the b ackground tempe rature (|Longaii) . 1 19661 : 
iPoolev fc Rvi3, [l968: Wall, 1990; Fc retti et al.l . l2004h . With 
the advent of new absolute measurements of the radio back- 
ground, coupled with radio source counts to ever increasing 
depths, the topic has un dergone something o f a revival. 

Recently a paper bv lGervasi et al.l (|2008l ) attempted to 
obtain fits to the source count data across a range of fre- 
quencies from v = 150 to 8440 MHz. From their fits, which 
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ranged from 1 /^Jy to 100 Jy, they were able to integrate the 
source counts to obtain an estimate of the sky brightness 
temperature contribution at each of the frequencies. They 
determined a power-law sky brightness temperature depen- 
dency on frequency with a spectral index of -2.7, which is in 
agreement with the frequency dependence of the fiux emit- 
ted by synchrotron dominated steep-spectrum radio sources. 
These estimates were used to interpret absolute measure- 
mon ts of the radio sk y brightness by the TRIS experiment 
( Za nnoni et all . l2008h . 

More recently the results of the 20 06 ARCADE 
2 balloon -borne experim ent were released (|Seiffert et all . 
[2009 ; Fi xsen et al.l . [2OO9I ). This instrument provided abso- 
lute measurements of the sky temperature at 3, 8, 10, 30, 
and 90 GHz. These results showed a measured temperature 
of the radio background about 5 times greater than that cur- 
rently determined from radio source counts, with the most 
notable excess of emission being detected at 3 GHz. Since 
most systematic effects explaining this emission were ruled 
out, we are left with the question of whether it could be 
caused by some previously unknown source of extragalactic 
emission. 

It was suggested by Seiffert et al. in the ARCADE 2 re- 
sults paper that this excess emission may be coming from the 
sub-/.tjy range. One might imagine an unknown population 
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of discrete sources existing below the flux lim it of current 
survey s. This issue was further examined by ISingal et al.l 
1)2010! ). Taking into account that a class of low flux sources 
must extend to ~ 10~^/^Jy (at 1.4 GHz), they concluded 
that this emission could primarily be coming from ordinary 
star-forming galaxies at 2: > 1 if the radio to far-infrared 
observed flux ratio increases with redshift. 

Before looking for radical causes of this emission, it is 
worth reexamining the observed radio source data to see if 
the ARCADE 2 result really does differ from what is ex- 
pected. To do this we derive new estimates of the source- 
integrated CRB at various frequencies and derive formal er- 
ror estimates for each. In Section [2] we describe the source 
count data used, together with our procedure and results 
for fitting the observed radio data. In Section [3] we present 
our estimates for the background sky temperature contribu- 
tions and the analysis of the uncertainties associated with 
these estimates. In Section|4]we compare our results to those 
obtained by the ARCADE 2 and TRIS collaborations. 



2 THE RADIO SOURCES AND THEIR 
COUNTS 

2.1 The Data Set 

Radio source counts at lower frequencies have been avail- 
able since the 1960s. There are many compilations of radio 
source counts available , particularly i n the l ast decade (e.g . 
Fomalo nt erall. 12 002": 'Bo ndi et all l2003l : iHopkins et"al] . 
2003.: Prandoni et al.. . ,2006i '). More recently deep continuum 
surveys at higher frequencies have become available, and 
with the use of newer technologies, have dramatically in- 
creased the amount and quality of data. The data used in 
this paper are from continuum surveys carried out from 1979 
to 2009 (see De Zotti et al. 2009, Sirothia et al. 2009). We 
used source count distributions from 150 MHz to 8400 MHz, 
with the individual frequencies covered being = 150 MHz, 
325 MHz, 408 MHz, 610 MHz, 1.4 GHz, 4.8 GHz and 8.4 GHz. 
References for all number counts used can be found in Ta- 
ble [1] 



2.2 The Number Counts Fit 

For fitting the source count data we opted to use a fifth or- 
der polynomial. A third order polyno mial was used in source 
count fitting bv lKatgert et al.l ()l988l l and a si xth order poly- 
nomia l fit to the 1.4 GHz d a ta wa s used by iHopk ins et al.l 
(|2003D , while ICervasi et al.l (|2008l ) used simple power-law 
fitting. Polynomial fits are simpler than some other choices 
of function, but still allow for fitting of different features 
in the data, such as the upturn at the low fiux end seen 
at some of the frequencies (where we note that an addi- 
tional sub-mjy peak could make a substantial contribution 
to the background). We chose a fifth order polynomial as it 
is high enough order to account for the features seen in the 
1.4 GHz data. Going to higher orders creates unneccesary ex- 
tra parameters while not improving the by a significant 
amount. Our empirical fits are performed on the Euclidean- 
normalized counts, i.e. F{S) = S^'^{dN/dS), with S being 
the flux density in Jy, using the polynomial with parameters 

F(S) = Ao + AiS + A^S^ + AaS'-" + AiS"- + AsS"". (1) 



The fitting is initially performed using a minimiza- 
tion routine. The minima are then used as starting 
points in a Monte Carl o Markov Chain, or MCMC approach 
(|Lewis fc Bridlel . l2002l ). which is used to refine the fits and 
obtain estimates of uncertainty. More details on the MCMC 
method can be found in section [321 The best fit values for 
all the parameters at each of the frequency bands can be 
found in Table [2] along with x^ values for each fit. The data 
and the best fit lines are plotted in Fig. [T] which shows the 
Euclidean normalized data, as well as the normalized re- 
sults. These (dN/dS) (surface brightness per logarithmic 
interval in flux density) plots are included to show where 
the peak contributions to the background arises. The right- 
hand panels in Fig. [T]show that the bulk of the background 
comes from relatively bright radio sources, with S ~ 1 Jy at 
the lowest frequencies to tens of mjy at the highest frequen- 
cies. But there is a significant, and still poorly characterized, 
contribution from much fainter sources. 

Table[2]shows that the x^ values of the fits are generally 
good, with all but one of the reduced x^ values being below 
2. The exception is for the 1.4 GHz data set, with a x^ of 
over 20 per degree of freedom. To obtain anything like a 
reasonable x^ ""^^ would have to increase the errors by a 
factor of four. It is worrisome that the 1.4 GHz compilation 
is the one with the most available data. As can be seen in 
the plot, there are many data points that are inconsistent 
with each other, even with the relatively large error bars. 

There are clearly systematic differences between differ- 
ent surveys at 1.4 GHz, particularly at the faint end. In the 
^Jy range it is difficult to obtain reliable counts, as this 
range is close to the natural confusion limit of mos t radio 
surveys (|Condon fc Mitchelll . 1 19841 : IWindhorst et al.1 . Il985l ) 
and hence the level of incompleteness may be incorrectly es- 
timated in some surveys. Moreover, at the bright end there 
are significant and systematic sources of error introduced 
when attempting to correct for source e xtension and surfac e 
brightness limitations (see discussion in lSingal et al.l . [20101 ). 
In addition to these effects, sampling variance (enhanced by 
source clustering) can lead to differences in counts for small 
fields. All of these systematic effects make it difficult to as- 
sess robustly the uncertainties in the derived CRB, as we 
discuss in the next section. 



3 CONTRIBUTION TO SKY BRIGHTNESS 
TEMPERATURE 

3.1 Integration of Radio Counts 

We integrate best-fit polynomials to obtain the contribu- 
tion from the sources to the sky brightness. To do this we 
integrate the function S{dN/dS) for each data set only in 
the range where data are available. We make this conserva- 
tive choice to avoid extrapolating at the very low and high 
flux density ends. Because of this our estimates of the sky 
brightness should be seen as lower limits. Thus to estimate 
the intensity we integrate 



(2) 



where Smin and Smax are different for each frequency. Once 
the intensity is determined we use the Rayleigh- Jeans ap- 
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Figure 1. Left: DifTercntial source counts Euclidean normalized and multiplied by c, with c = 1000,100, 10, 1, 1, 1, and 0.1 for u = 
150, 325, 408, 610, 1400, 4800, and 8400 MHz respectively. Right: 5^ normalization to show where the contribution to the sky brightness 
temperature is largest. The counts are a compilation from many different surveys, listed in Table [T] Solid lines are best fit polynomials. 



Table 1. References for the Extragalactic Radio Count Data Compilation 



Frequency 



References 



150 MHz 
325 MHz 
408 MHz 
610 MHz 

1.4 GHz 



4.8 GHz 
8.4 GHz 



Hales et al. 1988 : McGilchrist et al. 1990 . 

Owen fc Morrison 2008 : Port et al. 1988: Sirothia et al. 2009 . 
Benn et al. 1982 : Grueff 1988: Robertson 1973 . 
Bondi et al. 2007 : Garn et al. 2008 : Ibar et al. 2009 . 
Katgert 1979 : Moss et al. 2007 . 

Bondi et al. 2008 : Bridle et al. 197 2: Cilie gi et al. 1999 : Fomal ont et al. 2006 : Gruppioni et al. 1999 : 
Hopkins et al. 2003 : Ibar et al. 2009: Mille r et al. 2008 : Mitch ell fc Condon 1985 : Owen fc Morrison 2008 : 
Richards 200 0: Seym our et al. 200 8: Whit e et al. 1997 . 

Altschuler 1986 : Donnelly et al. 1987: Fo malont et al. 1984: Gregory et al. 1996 ; 
Kuehr et al. 1981: Paulinv-Toth et al. 1980 : Wrobel &z Krause 1990 . 
Fomalont et al. 2002 : Henkel fc Partridge 2005 : Windhorst et al. 1993 . 



proximation to convert it to a brightness temperature, 

r(i^) = /M£, (3) 

where k is the Boltzmann constant. The results from the 
integration at each of the seven frequencies are listed in Ta- 
ble O 

After obtaining these conservative estimates we next 
investigate the effect of reasonable extrapolations on the re- 



sults, with the limits of integration broadened to 10~^ Jy 
and 10^ Jy for Smin and S'max, respectively. The 1.4 GHz 
data set has the most extensive coverage across the flux 
density range. For this reason, we manually extrapolate the 
curve for the 1.4 GHz data set and integrate to get a new 
estimate for the background temperature. To extend the lim- 
its of integration for the 1.4 GHz data the end behaviour of 
the polynomial is constrained with the assumption that the 
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Table 2. values for best fits at each of the frequencies 



u Degrees of Aq A\ A2 A3 A4, A5 

MHz Freedom 



150 


68 


45 


6.58 


0.36 


-0.65 


-0.19 


0.26 


0.099 


325 


59 


34 


5.17 


0.029 


-0.11 


0.36 


0.17 


0.20 


408 


66 


44 


4.13 


0.13 


-0.34 


-0.003 


0.035 


0.01 


610 


75 


59 


3.02 


0.71 


0.97 


0.91 


0.28 


0.028 


1400 


4230 


196 


2.53 


-0.052 


-0.020 


0.051 


0.010 


-0.0013 


4800 


32 


47 


1.95 


-0.076 


-0.15 


0.020 


0.0029 


-0.00079 


8400 


41 


29 


0.79 


-0.10 


-0.23 


-0.051 


-0.019 


-0.0029 



Table 3. Values of the integrated sky brightness and tempera- 
ture contribution from radio source counts for different frequency 
bands. The uncertainties are la limits determined from Markov 
chain polynomial fits to the data. The high and low extrapolations 
are discussed in the text. 
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Figure 2. Integration results and best fit power-law from Equa- 
tion H] 



counts fall off beyond the low-intensity end of the data. Ar- 
tificial points are added in this region and their positions 
varied until a reasonable fit to the data is achieved. This 
procedure is repeated with both a steep and shallow roll-off, 
to obtain high and low background estimates. These slopes 
are chosen to be the most reasonable steep and shallow es- 
timates, with the x^s being a factor of 5 and 7 greater than 
the best fit to the data alone. The best fits for the extrap- 
olations can be seen in Fig. [31 The higher estimate could 
have been allowed to have an even shallower slope, therefore 
allowing for an even higher background estimate; however 
anything much shallower than the chosen fit would have 
values several times larger again. This fact makes any shal- 
lower fits an unreasonable choice. The steep slope estimate 
for the 1.4 GHz data ends up giving nearly the same result 
for the background temperature as the unextrapolated es- 
timate. This is because the unextrapolated estimate has a 
rising low-intensity tail (Fig. 1), while the extrapolated es- 
timate (Fig. 3) has limits of integration extended with the 
roll-off procedure, whose end behaviour is controlled to pro- 
duce a steep downturn beyond the available data. 

From our conservative estimate of the 1.4 GHz back- 
ground a power-law is fit to the temperatures. This takes 
the form of 

^(^^-^(liSh;)'' W 

where A is the power-law amplitude, and /3 the index. We 
set A to the 1.4 GHz value of 0.110 K, while Monte Carlo 
Markov Chains (section l3.2|) were used to find the best value 



of /3 = -2.28 ± 0.1. The results of this power-law fit can be 
seen in Fig. (2] The fit is high for the two highest frequency 
points, primarily as a result of the 1.4 GHz data, which has 
the most flux density coverage and most available data. It 
should be noted that this is just a phenomenological fit; we 
see little merit in adopting a more complex model such as 
a broken power-law just to satisfy two extreme data points, 
particularly when this power-law is only used to obtain ex- 
trapolation estimates. 

With the high and low extrapolation estimates from the 
1.4 GHz data, we use Equation|4]to obtain estimates for the 
other frequencies. The results of the extrapolated estimates 
are given in Table [3] As even these reasonable extrapolations 
can change the background estimates by about a factor of 
two, it is clearly important to push counts at all frequencies 
to fainter levels. 



3.2 Uncertainty — Monte Carlo Markov Chains 

To investigate the uncertainties thoroughly, we carry out our 
fits with Monte Carlo Markov chains for eac h of the data 
sets, using CosmoMC (|Lewis fc Bridld , I2OO2D as a generic 
MCMC sampler. The function is sampled for each set 
using the polynomial in equation 1, which is then fed to the 
sampler to locate the minimum. Each of the six parame- 
ters of the polynomials are varied for each step of the chain 
and the chains are run with 500,000 steps. CosmoMC gener- 
ates statistics for the chains, including the minimum , the 
best fit values for each of the parameters, and their uncer- 
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Figure 4. 100 Markov Chain polynomial fits generated for the 
1.4 GHz data set. Greyscale indicates relative probability, with 
the solid black line being the best fit curve. 



3.3 Comparison with Previous Estimates 

Over the years there have not been many estimates 
of the CRB made using source count data |Longaii , 



Il96 6l: iPoolev fc Rvlel Il968| - IWalll . Il990l : iFeretti et all . 12004 
Gcrvasi et al.l . l2008l ). And even within this small list, the fre- 
quencies covered were rather limited and uncertainties not 
always quoted. It is important to see how our e stima tes com- 
pare with these previous estimates. Longaid (|l966h gives a 
value for Tirs = 23 ± 5 K. IWalll (|l990l ) lists estimates of 
T408 = 2.6 K, T1.4 = 0.09 K, and T2.5 = 0.02 K. Our resuhs 
are in agreement with these earlier estimates to within ± 
2(7. The values for source contributions from Gervasi et al. 
(2009) tend to be a little higher than ours, the differences 
being traceable to choices made for the limits of integration 
and for the parameterized form for the fits. 

The ARCADE 2 experiment reported an excess of emis- 
sion beyond what we and others have estimated from source 
counts, with the excess largest at 3.4 GHz. We have also 
considered much lower frequencies in this paper than the 
3.2 GHz detection limit of ARCADE 2. However, it is possi- 
ble to calculate what temperatures would be expected using 
the best fit to the ARCADE 2 data: 



T{u) = To + A\ 



(5) 



VlGHz/ 

Here Tb is the CMB base temperature, and the best fit 
values for the para meters are P = -2.56 and A = 1.06 
(|Seiffert et al.l . |2009| ). Measurements from the TRIS exper- 
iment were performed at v = 0.6, 0.82, and 2.5 GHz, and 
compared with the Gervasi et al. source contribution calcu- 
lations are within 3% at 0.6 GHz and 50% at 2.5 GHz 

The quantities detected by or extrapolated from AR- 
CADE 2, those estimated from counts by Gervasi et al. 
(2008), the measurements from the TRIS experiment, as 
well as our current estimates are shown in Fig. (6] Here it can 
be seen that the ARCADE 2 absolute measurements lie far 
above both source estimates and TRIS measurements, par- 
ticularly at lower frequencies. Clearly, the excess detected 
around 3 GHz would correspond to a large excess at lower 
frequencies if the power-law continued. 



tainties. As an example, Fig. [4] shows different polynomial 
fits tested by the MCMC and their relative probability for 
the 1.4 GHz data set. 

Histograms of the chain values for the background tem- 
perature are shown in Fig. [S] From the width of these his- 
tograms we are able to measure the uncertainty in our es- 
timates for the background temperature, taken here as the 
68 percent area values, fully accounting for the correlations 
among the parameters in the polynomial fits. The la uncer- 
tainties are listed in Table [3] 

Most of the histograms are fairly Gaussian, which is a 
refiection of the quality of the data. Frequencies with good 
data around the peak contribution (in the right-hand panels 
of Fig. 1) tend to have well-constrained background temper- 
ature values, e.g. at 408 MHz. However, there is a notice- 
able irregularity with the 325 MHz histogram. Because of 
the limited data available at 325 MHz, and with the peak 
area of contribution having little to no data, the histogram 
at this frequency does not have a well defined shape; the 
uncertainty is far from Gaussian. 



3.4 Systematic Errors 

We have considered several possibilities for systematic er- 
rors in exploring whether our results might be compatible 
with those from the ARCADE 2 experiment. The first of 
these is possible bias from source clustering. This can be 
an issue when dealing with surveys covering small areas, 
where one might get more field-to-field variations than ex- 
pected from Poisson errors. The two-point angular corre- 
lation function for NVSS and FIRST source s fits a power- 
law shape for separa t ions u p to at least 4° (|Blake fc Wall 
l2002l : [O verzier et al.l . l2003h . From this angular correlation 
one can estimate the fr actional variance of the counts 
JSeldner fc Feeble j. Il980l l. This procedure was carried out 
bv lde Zotti et al.l (|201o[ )" and has been taken into account in 
the errors provided and used in our estimations. 

Another effect that could influence our results is the fact 
that in some of the surveys used in our compilation the mea- 
surement frequency was slightly different from the nominal 
one, i.e. 5 GHz rather than 4.8 GHz. In such cases we scaled 
the original measurements to the nominal frequency using 
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Figure 5. Histograms from the Markov chains at each frequency. The background temperature was computed at each step in the chain 
and binned. At most frequencies there is a well defined value of the CRB with an approximately Gaussian distribution, while this is less 
true at 325 MHz in particular. 
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Figure 6. Integrated sky brightness temperature at each fre- 
quency from the estimates in this paper (open squares), Gervasi 
et al (2008, stars), TRIS measurements (Zannoni et al. 2008, di- 
amonds), and ARCADE 2 measurements (Seiffert et al. 2009, 
triangles). The dashed line at t he bo ttom represents the CMB 
temperature at 2.726 K l|Fixsenl . liooil ') . 

the assumed dependence of the source flux S{v) ~ v^^'"^ . 
This correction resuhs in neghgible change in the derived 
fits. 

An additional effect that could account for the differ- 
ence in the background temperatures is the possibility that 
some surveys have somehow missed extended high-frequency 
emission blobs which could integrate up to the required 
amounts. This seems an unlikely option, as such structures 
would have to be on degree scales or larger to escape detec- 
tion, and because if these structures have features above a 
certain brightness temperature then they would have been 
seen. 

Other possible effects to take into consideration for the 
uncertainties include: 

1. Calibration variations for different radio telescopes; 

2. Inaccurate determination of completeness corrections 
at the faint end; 

3. Contribution from diffuse emission from the Intergalac- 
tic Medium (IGM), Intercluster Medium (ICM), and the 
Warm-Hot Intergalactic Medium (WHIM) ; 

4. Missing low surface brightness emission from extended 
objects that are either large, or sources with extended com- 



ponents, or sources that are not detected if source surface 
brightness extends to low values. 

ISingal et all (|2010l ) provide a detailed discussion of 
items 3 and 4 as well as several other possibilities such as 
radio supernovae that could contribute to the CRB. We sus- 
pect that the most important effects are the first two items, 
particularly completeness at the faint end. 



4 INVESTIGATION OF A FAINT 'BUMP' IN 
THE COUNTS 

At 1.4 GHz, where we have the most data, our estimated 
background temperature plus a CMB baseline is 2.83 K ± 
0.02 K, while an extrapolation of the ARCADE 2 result gives 
3.17 K ± 0.01 K (error estimate from their measurement at 
3.2 GHz). This corresponds to a difference that is nearly llo 
away from our estimate. It has been suggested that this 
could be explained through an extra population of faint ra- 
dio galaxies, corresponding to a 'bump' in the Euclidean- 
normalized counts at flux densities near or below where th e 
current data are petering out (see also ISingal et al.l I2OIOI ) . 
We want to investigate how big this bump would need to be 
in order to explain the excess emission. 

We carried out two separate approaches for modelling 
such a bump, the outcomes of which can be seen in Fig. [T] 
Our first approach is a simple extension of the current counts 
with an upward trend below 10 /iJy, but one not quite as 
steep as the best fit line. To do this we simply added artificial 
data points past the lower flux density limit of the rest of the 
data in order to control the end behaviour of the fit line. We 
then investigated what was required to match the ARCADE 
2 results. The solid line of Fig. [7] shows the results of this 
fitting. It is the best fit to the data, allowing for a moderate 
upward slope in the faint end. When integrated from 10~® Jy 
the result is enough to account for the temperature reported 
by ARCADE 2. 

Our second method involved choosing a simple parabola 
with fixed width of a decade in log-S and variable position 
for the peak, and running a Markov chain that fit the height 
parameter that would integrate to give the amount of excess 
emission needed to match the ARCADE 2 result. We found 
that the peak of the bump could be at flux densities as high 
as 8.0 /iJy. 

It is relatively easy to produce a bump big enough to 
account for the extra emission while still fitting the rest of 
the data reasonably well, with either method. However, we 
do know that any such bump is constrained by the observed 
IR background, through th e IR-Radio correlation (see e.g. 
iHaarsma fc Partridgel Il998l ) . This correlation will have to 
be taken into account in any modelling of this faint fiux 
density bump so as not to overproduce the IR background. 
This essentially requires any faint radio population to be 
quite IR faint compared with known galaxy types. 



5 CONCLUSIONS 

We used source count data from v — 150 MHz, 325 MHz, 
408 MHz, 610 MHz, 1.4 GHz, 4.8 GHz, and 8.4 GHz to evalu- 
ate the contribution from sources to the diffuse cosmic radio 
background. Polynomials were fit to the data and integrated 
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Figure 7. The 1.4 GHz source count data. The thin soUd Unc 
gives the best fit to the data while having a moderate slope at 
the faint end and integrating to the necessary excess reported by 
ARCADE 2. The thick solid Une is our best fit to the 1.4 GHz 
data from Section 2. The other three lines are bumps peaking at 
7.9 (dotted), 5.0 (dashed), and 3.1 (dot dash) /iJy, which produce 
the background temperature necessary to match the ARCADE 2 
results. On this plot the height of such a bump is proportional to 
ol/2 

peak 



to obtain lower bound estimates at each frequency for the 
sky brightness temperature. In addition, we also extrapo- 
lated our fits beyond the limits where data are available 
using reasonable assumptions for how the curves behave in 
those regions. We then used Monte Carlo Markov Chains to 
obtain estimates of the uncertainties of the temperature es- 
timates at each frequency and also considered other possible 
sources of uncertainties that could affect the results. 

Our estimates are considerably lower than the measure- 
ments of ARCADE 2, even when taking into account the un- 
certainties or extrapolations. We considered the possibility 
that the excess emission comes from a bump in the source 
counts in the /xJy range at 1.4 GHz. We used modelling to 
see how large such a bump must be in order to obtain the 
necessary contribution to the background. We saw that a 
bump could exist in this range, peaking at fluxes as bright 
as 8 ^Jy, and could integrate up to the excess emission of 
± 320 mK, with a height that is consistent with the data. 

We still have no direct evidence that such a new popu- 
lation exists, and so further investigation into the faint end 
of the counts is needed. The infrared and radio connection 
could be used to test this idea through use of signal stacking 
and by examining different possible luminosity functions to 
look at the evolution of such a population. The final answer 
may only be reached when source count data become avail- 
able in the /ijy range, perhaps in the era of the EVLA and 
eventually the SKA. 
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